Intensity correlations in resonance nonlinear magneto-optical rotation 
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We have studied the intensity correlations between two orthogonally linearly polarized compo- 
nents of a laser field propagating through a resonant atomic medium. These experiments have been 
performed in a Rubidium atomic vapor. We observe that the correlations between the orthogonally 
polarized components of the laser beam are maximal in the absence of a magnetic field. The mag- 
nitude of the correlations depends on the applied magnetic field, and the magnitude first decreases 
and then increases with increasing magnetic field. Minimal correlations and maximal rotation angles 
are observed at the same magnetic fields. The width of the correlation function is directly propor- 
tional to the excited state lifetime and inversely proportional to the Rabi frequency of laser field. 
These results can be useful for improving optical magnetometers and for optical field or atomic spin 
squeezing. 
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I. INTRODUCTION 

Quantum coherence effects have been a focus of 
research activities for the last two decades, because 
they may drastically change the optical properties of 
a medium. For example, electromagnetically induced 
transparency (EIT) [U 0K predicted and observed in CW 
and pulsed regimes d [j] , practically allows absorption 
to vanish. The nonlinear response of a resonant atomic 
medium at moderate optical intensities can be strongly 
enhanced by creating maximal coherence between long- 
lived sublevels of the ground state @, S H, 0, H- The 
corresponding steep dispersion results in the ultra-slow or 
fast propagation of light pulses [9|, lid 1 1 XI ] , which can pro- 
duce huge optical delays [T2| and can be used for drastic 
modification of the phase-matching conditions for Bril- 
louin scattering [l3j| . and four- wave mixing [bj ]. ft is 
possible to manipulate a coherent medium and produce 
optical pulses at rates faster than the relaxation rates of 
the medium [l5j], in order to develop bri ght sources for 
efficient generation of IR and FIR pulses [16} . Also, the 
nonlinear properties of such media are enhanced, allow- 
ing the implementation of quantum li ght storage [It], EH , 
nonlinear optics at a few photon level [l9( , and other non- 
linear effects |20l | . 

It is worth noting here that media with quantum coher- 
ence might be used for the sensitive measurement of mag- 
netic fields [2lj . Recent interest in nonlinear magneto- 
optical effects (NMOE), such as nonlinear magneto- 
optical rotation (NMOR) [5j or nonlinear polarization 
self- rotation (NPSR) [22|], has been fueled by possible 
applications to precision magnetometry [23l . [24 . l25l . I2H ] . 

The amplitude noise of optical fields increases after the 
fields interact with an atomic media (atomic excess noise) 
[H HE El • Several processes, such as the conver- 
sion of phase noise to amplitude noise and the four-wave 



mixing process, can be associated with the generation of 
atomic noise [11 HI, H . The strong coupling 
of optical fields in A-type atoms provides a way of control- 
ling amplitude fluctuations. The amplitude correlations 
have been studied in EIT 0,0 011 13 and NMOE 
[38| experiments. Observation of sub-Poisson statistics of 
the amplitude noise is reported in [34j . Amplitude corre- 
lations and anti-correlations are observed in two orthog- 
onally circularly polarized optical beams from the same 
laser (3(| and in orthogonally linear ly p olarized optical 
beams from two independent lasers (35[ . Power spectra 
of intensity fluctuations was studied in [37j ■ Noise spec- 
troscopy of nonlinear magneto-optical resonances in Rb 
vapor at a detection frequency 2.5 MHz with bandwidth 
30 kHz is performed in [381 ] , where the photon noise limit 
has been reached. 

Coherent effects may be used to reduce the noise level 
below the limit given by photon fluctuations. Several 
schemes for squeezing of the fields are proposed and stud- 
ied in EIT US M or NPSR HI S3, El El con- 
figurations. It has been demonstrated that the atomic 
excess noise can reduce squeezing in the fields [45j, and 
the addition of excess noise to quantum state of radia- 
tion is important and should be taken into account in 
experiments involving coherently prepared atoms. For 
instance, recently a considerable influence of the excess 
noise on light storage has been observed jig} . 

In this paper, we have considered the atomic noise cor- 
relations in an NMOR experiment with a rubidium va- 
por. In a previous paper [3t3 |. we studied the fluctua- 
tions between left and right circularly polarized beams, 
which are the normal modes of the system. Here, we 
have studied the fluctuations between two orthogonally 
linearly polarized beams, which are not the normal modes 
of the system. Our current experimental configuration is 
widely used in optical magnetometry measurements 
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FIG. 1: Experimental setup (a) and energy levels (b). The 
ground state Zeeman splitting of Rb 87 is 0.7 MHz/G. 



and therefore furnishes an ideal system for studying these 
fluctuations. We have observed the variation of the corre- 
lation function G^ 2 '{t) (t is the delay time) as a function 
of the magnitude of an applied magnetic field. 

The paper is organized as follows. In Section II, we 
describe the experimental setup and the results. In Sec- 
tion III, we present a theoretical model, which gives us 
an explanation of the results. Finally, in Section IV, we 
present our conclusions. 



II. EXPERIMENTAL SETUP AND OBTAINED 
RESULTS 

The schematic of our experimental setup is shown 
in Fig. la. An external cavity diode laser (4?| (the 
linewidth is less than 1 MHz) is tuned to the center of the 
Doppler broadened D\ line (transition bS\/ 2 (F = 2) — > 
5P\/2{F = 1)) of 87 Rb (see Fig. lb). After a polarizing 
beam splitter (PBS), the laser beam (of diameter 0.1 cm 
and power 0.9 mW) enters a glass cell (of length 7.5 cm) 
filled with an atomic Rubidium vapor (with natural iso- 
tope abundance and atomic density 10 12 cm~ 3 ). 

The cell is placed inside a two-layer magnetic shield 
that protects the cell from uncontrolled environmental 
magnetic fields. A solenoid is installed inside the shield to 
create a longitudinal magnetic field B. A half-wave plate 
is placed between the two PBSs to rotate the polarization 
of the output beam to 45° with respect to transmission 
axis of the second PBS. This angle is set by making the 
two beams have equal intensities at a large one-photon 
detuning from the atomic resonance in the cell at room 
temperature. Two orthogonally linearly polarized beams 
from the second PBS are focused on identical fast photo- 
detectors (PDs) with a frequency bandwidth of 75 kHz- 



1.2 GHz. The optical path lengths for both beams and 
length of cables between the PDs and the data recorder 
are chosen to be the same to avoid a time delay between 
the signals in two registration channels. The signals from 
PDs are recorded by a digital storage oscilloscope (DSO) 
with a bandwidth of 100 MHz. 

We have also studied the transmission and the polar- 
ization rotation of the beams as a function of the two- 
photon detuning by scanning the magnitude of the lon- 
gitudinal magnetic field B at a frequency of several Hz. 
For this, the fast photo-detectors have been replaced by 
low-frequency detectors to monitor the output from the 
second beam splitter. Then, using the recorded signals, 
Si and 5*2, we calculate the transmission by 

(sbi + S02) ^ ^ 

and the polarization rotation by 

= arcsin (frf ) - (2) 

where Sqi and S02 are the signals from the photo- 
detectors without atomic resonance absorption. 
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FIG. 2: Magnetic field dependence of the EIT (Fig. 2a), po- 
larization rotation (Fig. 2b) and correlation function G' 2 '(0)) 
(Fig. 2b). In Fig. 2(c) the circles represent experimental data 
and the solid curve is a visual guide for eyes. 



We start the presentation of the results with a demon- 
stration of the narrow EIT resonance. The narrow EIT 
resonance is related to the "dark state" of Rb atoms 
formed by the action of two laser fields. The transmis- 
sion T versus magnetic field B is shown in Fig. 2a. In 
a magnetic field, the linear polarization of the laser field 
rotates, and the polarization rotation dependence on the 
external magnetic field is shown in Fig. 2b. The width 
of nonlinear magnetic optical resonances is determined 
by the power broadening of the two-photon transition, 
which for our conditions is narrower than the natural op- 
tical width (6 MHz). The observed background in Fig. 2b 
is due to the linear rotation of polarization. 

To study the fluctuations of an optical field transmit- 
ted through a dense Rb vapor, we have registered the 
time dependent signal fluctuations 551,2 (i) of the- two 
optical beams after the second PBS. We record these sig- 
nal fluctuations for the two output beams in a 10 [is 
time-window for different magnetic fields, and then we 
calculate the normalized correlation functions given by 



(&gl(t)6g 2 (t + T)) 



(3) 



where SSi^(t) are the time dependent fluctuations of the 
two beams, and stochastic averaging [5l[ denoted by an- 
gular brackets is defined as (Q(t)) = 1/T f* +T Q(t')dt'. 

We are unable to detect any intensity fluctuations 
when the cell has been removed from the laser beams. 
Of course, it is well-known that diode laser radiation 
possesses low intensity noise. In our experiments, the 
photo-detector noise dominates, as in Refs. [36l l37l] . 

In the presence of a resonant medium, however, the 
situation changes, and the phase noise of a diode laser 
is transformed into intensity fluctuations. The magnetic 
dependence of the correlation function at zero time delay, 
G^ (0), is shown in Fig. 2c. As one can see, a correlation 
magnitude of close to 0.9 is obtained at zero magnetic 
field. 

Increasing the magnetic field up to 0.2 G results in 
decreasing the correlation function magnitude to ~ 0.2. 
But further increase of the magnetic field leads to a re- 
vival of the correlation, and it reaches ~ 0.7 at 0.8 G 
field strength. We underline here that this behavior is 
different from that observed in [3(|. The key difference 
is that, in this paper, we study the fluctuations of the or- 
thogonal linear polarization components, which are not 
normal modes because of the Faraday effect. 

By examining the polarization angle rotation depen- 
dence shown in Fig. 2b, one can conclude that the maxi- 
mum correlation is obtained when the magnetic field and 
the polarization rotation in the atomic medium are close 
to zero. The correlation functions G^ (r) for three dif- 
ferent magnetic fields are shown in Fig. 3. 

The signals from the photo-detectors in time intervals 
of 200 ns are shown in the inset boxes in Fig. 3a,b,c. The 
signals are proportional to the laser beam power with a 
slope of 500 V/W. The vertical range of the plot is 4 mV. 
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FIG. 3: Intensity correlation function G' 2 ' (t): curve (a) 
recorded at zero magnetic field, curve (b) at a magnetic field 
of 0.18 G and curve (c) at a magnetic field of 0.9 G. Also 
signals from photo-detectors are shown in boxes. The wave- 
forms are recorded in time intervals of 200 ns and amplitude 
intervals of 4 mW. 



The amplitude variations of the signals are practically the 
same as in Fig. 3. 

The temporal behavior of the signal is modified by the 
applied magnetic field. The signals shown in Fig. 3a and 
3c are correlated, and the ones shown in Fig. 3b (box) 
are not. We observe that at low and high magnetic field, 
fluctuations are strongly correlated. At an intermediate 
magnetic fields, the peaks of correlation function have 
less magnitude. The correlation function behaves as if it 
changes sign around zero time delay. Possibly different 
frequency components of the signals can be correlated or 
anti-correlated at these magnetic fields. For the depen- 
dence shown in Fig. 2c, we have selected the magnitudes 
of the correlation function at zero time delay. 

We note that the inverse width of the peaks is the order 
of 27r -20 MHz, and it is comparable with the width of the 
saturation resonance [H, |49|, , and it is much broader 
than the ground state relaxation rate. The inverse width 
depends on the excited state decay rate 27r • 6 MHz and 
the optical excitation rate. The influence of one-photon 
optical saturation on MNOR is discussed in [j| . 
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FIG. 4: A simplified three-level scheme. Splitting of the 
ground state is u> c t = 2q^lbB, where g is the Lande factor, 
Hb is Borh's magneton, and B is the magnetic field. 



III. THEORETICAL APPROACH AND 
DISCUSSION 

The following is a simplified theoretical description of 
our system treating the laser fields classically and using a 
density matrix for the atomic response. The laser beams 
are in resonance with a three-level medium as depicted 
in Fig. 4. 

The Hamiltonian of the atom is given by 



h = nn^\a)(b\ + nn + \a)(c\ + h.c. 



where fl± = p c ,bE±/Ti are the Rabi frequencies of lcft- 
and right-circularly polarized beams; p c j> and E± are the 
corresponding dipolc moments of the atomic transitions 
and the electric fields. The density matrix equation is 
given by 



1 



9 : 



[H,p]--(Tp + pT), 



(5) 



that the left- and right-polarized beams are the normal 
modes of the current system, and their polarizations do 
not change while they are propagating through the cell. 

The intensities of the optical beams propagating 
through the cell fluctuate. The corresponding correla- 
tion function is defined as 



G< 2 )(r) 



(SI-(t)SI+(t + T)) 



(7) 



y/([6I-(t)]*)([6I + (t + T)]*) 

Intensity fluctuations have been experimentally and the- 
oretically studied in [H, [3(| , and they are related to the 
atomic responses 



51- 



%p ab n*_, 51+ ~ Qp ac fl* 



(8) 



The corresponding atomic coherences in the three-level 
system can be found from the solution of the density ma- 
trix equations. Assuming the phase diffusion of optical 
fields is a slow process, the atomic coherences are given 
by 



Pab = -l- 



.n ba fl- + p c .b^+ 



+ p cfe fL 



FcbPcb = iPcaQ- - iPabft+, 



( 4 ) where r afc = j ab + i{uj ab — V-)\ T ca = j ca - i{u a 



r, 



b = leb + i{LO c b — V-+V^ 



la/3 



(9) 



(10) 



us a p are the atomic frequencies; a and j3 are labels for 
atomic levels a,b,c; V—(t) = v+(t) are the instantaneous 
frequencies of laser radiation in both beams having or- 
thogonal polarizations. Assuming f2_ = Q + = fl, and 
defining 5 = uj ab — i>-— uj cb /2, we solve Eq. IjlOp with re- 
spect to p c b- Then, substituting the solution into Eqs. © 
gives us the following: 



where T is the relaxation matrix, introduced to describe 
relaxation processes in the atomic medium [51]. The 
equations for field propagation are 



dz 



= -iVbPab, 



<9Q 4 



= -IVcPa 



(6) 



where r\ b — V-Np\j \2%€qc), r\ c — v + Np^,/{2Ktoc) are 
the coupling constants, v± are the frequencies of circular 
polarized fields, p b ,c arc the dipolc moments of the corre- 
sponding transitions, N is the density of medium. Note 



_ . r cfc (r - iS)n ab + n cb Q 2 

Pah ~ 1 r cfc (<P + r 2 ) + 2rf> 2 ' 



T cb (T + iS)n ca + n cb tt 2 

1 r cb (5 2 + r 2 ) + 2rr> 2 



ft. 



(ii) 



(12) 



The set of equations for populations can be obtained by 
substituting Eqs. (|11I12[) into Eq.© for the appropriate 
atomic populations: 



/ r cfe (r - iS) T bc (T*+iS) \ 2 
JaUa + ~ 1 ~ \\l n ab 



cb 



cb 



J- + J- ) Vl A n cb = 0, 

r c6 r; 



T cb (T + i6) T bc (T* - iS) \ 2 
j a n a + ^ H ^- II n a 



cb 



1 



cb. 



i 



+ s- + n n bc = o, 



r 



(13) 
(14) 



cb , 



5 



where T cb = T cb (S 2 +T 2 ) + 2TQ 2 . By introducing A b = 
A - SB, A c = A + SB, and 

/r^ + r^r*\ 2 B= M 
V r* cb J \r cb r*J 

(15) 



(16) 



the set equations for the populations can be re-written 
in the compact form 




j a n a + A b n ab + Cn cb = 0, 



7 a n a + A c n ac + Cn bc = 0, 



(17) 



(18) 



and solved (see Appendix |A|) . The solution for the pop- 
ulation difference in level b and c is given by 



n cb 



2-faBS 



3A b A c + C(3A b + 3A C + 4 7q ) + j a (A b + A c ) ' 



Simplifying Eqs. (|15|16|) . we obtain A ~ -f cb , B 



C~ and 

7 ' 



n cb 



AS 

w 



(19) 



(20) 



The intensity fluctuations are determined by atomic co- 
herences, which are 



7 c6 7 + AS 

SPab — fFTi ' ^Pac 



7 |n| 



7 cb 7 - AS 



(21) 



where the field phase fluctuation is related to the fre- 
quency deviation S that is restricted by the EIT win- 

q2 

dow [36], S . 

7 

We can clearly see from Eq. (|2~lj) that the intensity fluc- 
tuations have two contributions: the first originates from 
absorption, which is the same for both modes, and the 
second originates from the Raman term appearing from 
the population difference n cb , which has opposite signs 
for these two modes. One field is amplified due to absorp- 
tion of the second field. Depending on the one- and two- 
photon detuning, A, the intensity fluctuations occur in 
phase if ^ cb ^ 2 /Vi 2 » Aor out of phase if 7 c& 7 2 /f2 2 <C A. 
The last condition gives rise to the anti-correlations. 

Let us note here that at the EIT condition, A = 0, cor- 
related intensities can be also considered from the point 
of view of matched pulses [52| . The details of a study of 
switching between correlations and anti-correlations can 
be found in [33|, [36| 

Now we are ready to consider intensity fluctuations 
between the orthogonally linearly polarized beams in the 
current experiment. Linearly polarized light can be re- 
garded as a linear combination of left- and right-circular 



polarized light. The light with left- and right-circular 
polarizations do not change their state of polarization 
while propagating through the gas cell. We previously 
studied the correlation of intention fluctuations between 
the normal modes [36| . 

In the current experiment, to detect rotation of po- 
larization, we use the following scheme. After the po- 
larizer, the beam propagates through the cell, and then 
after a half- wave plate, the polarization plane rotates 45° 
degrees and is split by a PBS into two orthogonally po- 
larized optical beams. The intensity of each beam is de- 
tected. Thus, the measured signals, Si and Si, are pro- 
portional to the intensities of the propagated circularly 
polarized components, which are given by 

1 



5i, 2 ~/i,2 = -(/++/- ± 2 sin 0), (22) 



and the polarization rotation is given by 



h-h 



h+h 



(23) 



Then correlation between S\ and S2 can be calculated 

by 

G ( 2 ) (T)= (ShWiit + T)) 

V([«iW] 2 >([W 2 (« + r)P) 

As is clearly seen from Eq. ([22"l) and Eq. ((24|) . the cir- 
cularly polarized beams are equally split into two linear 
polarizations. Thus, the two intensities of the linearly 
polarized components are correlated. The only contribu- 
tion that decreases the correlations is due to the terms 
depending on the rotational angle, which have different 
signs for orthogonal polarizations in Eq. (|2"2"]) . 

The intensity fluctuations should be small for the com- 
ponents that are not perfectly anti-correlated. But as 
we have shown in a previous paper [36], increasing the 
magnetic field causes these components to become anti- 
correlated. On another hand, increasing the magnetic 
field increases the rotation at first, but then once the 
splitting becomes bigger than the EIT width, the rota- 
tion angle decreases, restoring the correlations between 
the intensities of the linear components. 

One can see from Eq. (pM)) that if the rotational angle in 
a magnetic field is small, the beams mainly consist of the 
sum of intensities for left- and right-circularly polarized 
beams, and thus the fluctuations are correlated. Anti- 
correlations come from the terms having different signs 
due to Faraday rotation in the magnetic field. So, at zero 
magnetic field, we have only correlations. Note here that 
the sum of intensities is a major contribution and always 
contributes to correlations between intensities. But for 
stronger magnetic fields, due to the decreasing of the an- 
gle of rotation, this term vanishes, and only correlation 
between beams is observed. 

In the experiment, on the other hand, we observe the 
decreasing of correlation with increasing magnetic field 



intensity. It looks similar to the results obtained for cir- 
cular polarization, but the physics is completely different. 
Actually, the strong decrease of correlations observed in 
the experiment occurs because of strong anti-correlations 
for higher two-photon detunings. 

For larger two-photon detuning, the nonlinear 
magneto-optical rotation of polarization becomes smaller 
and intensity correlations restore their correlations. De- 
fine 

I+ = I + i+, I- = I + i-, (25) 

where i± are the intensity fluctuations of the beams I±, 
correspondingly. Then, using y/1 + z ~ 1 + | — 
introducing x = i+ + i- and s = i + — i_ (note that 

M- - T = 4 *+*- _ (*+ + l -f = -(*+ - i-f = -s 2 )' 
we can rewrite Eqs. (|2"21 as 

s 2 

h,2 = 2I +x± (2I Q + x- — ) sin (26) 

4io 

I 



In the last equation, one can see that for small mag- 
netic fields, when the rotation angle is small, and the 
beam intensities are correlated [36j . and for strong mag- 
netic fields, when the rotation angle is also small and the 
beam intensities are anti-correlated [36j, the correlation 
function equals to unity. These results are independent 
of correlations or anti-correlations between the beams. 
Correspondingly, in the intermediate case, we have cor- 
relations that are less perfect. 

The foregoing discussion furnishes a qualitative de- 
scription of the correlation behavior we have observed 
in this experiment. To obtain quantitative agreement, 
we would need to take into account the intensity fluctua- 
tions more accurately, beyond the small variations as we 
have assumed here. We will present the results of such 
an analysis elsewhere. 

IV. CONCLUSIONS 

We have studied atomic noise correlations in a nonlin- 
ear magneto-optical rotation experiment with Rubidium 



6 

Thus we obtain 



6h, 2 =h,2-{h,2)=x{l± sin 0) T ./ ' sin (27) 



and 



((h-(h)) 2 ) = (x 2 )(l+sin0) 2 + ^ } J* } sin 2 (28) 



Finally, we can calculate correlation function G^ 2 \ de- 
fined by Eq. flHJ) , 



(29) 



I 

atomic vapor by using broadband detection. The corre- 
lations between the orthogonally polarized components 
of the laser beam are maximal in the absence of a mag- 
netic field. The width of the correlation function peak is 
proportional to the excited state lifetime and the inverse 
Rabi frequency. When a longitudinal magnetic field is 
applied, the correlations first decrease and then increase. 
The minimal correlations and the maximal rotation an- 
gles are observed at the same magnetic fields. These re- 
sults can be useful for improving optical magnetometers 
and for squeezing optical fields or atomic spins. 
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Mx 2 
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APPENDIX A: SOLUTION OF SET OF THE „ _ A b A c + C{A b + A c + 2 7a ) + A b ^ a 



DENSITY MATRIX EQUATIONS 3A b A c + C(3A b + 3A C + 4 7a ) + ^ a (A b + A c ) 

(A5) 

The set of equations from populations Ea. (|T?|T8j) is and population differences can be presented in the fol- 
given by lowing form 

j a n a + A b n ab + Cn ch = (Al) 



la {A c + 2C) 

n a b ■ 



7a n Q + A c n ac + Cn bc = (A2) 

where A b , A c , C are defined above (see Eqs. (|f 511 6[) ). 
Then solution of Eqs. (|AllA2|) are the following 



3A b A c + C(ZA b + 3A C + 4 7o ) + la {A b + A c ) ' 

(A6) 



la {A b + 2C) 



A b A c + C(A b + A c ) ZA b A c + C{ZA b + ZA c + A la ) + la {A b + A c y 

Ua ~ 3A b A c + C{ZA b + 3A C + 4 7a ) + la {A b + A c ) ( A7 ) 

(A3) 

A b A c + C{A b + A c + 2 7a ) + A cla la {A b ~A c ) 
n b = 777—, — ; a — r^n — r~A — r - ; — n — r^rr n cb ~ 



3A b A c + C{3A b + 3A C + 4 7a ) + la {A b + A c ) 3A b A c + C(3A b + 3A C + 4 7a ) + la {A b + A c ) ' 

(A4) (A8) 



